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Abstract 



The behaviour of Jacobi fields along a time-like geodesic running 
into an isotropic singularity is studied. It is shown that the Jacobi 
fields are crushed to zero length at a rate which is the same in every 
O , direction orthogonal to the geodesic. We show by means of a counter- 

example that this crushing effect depends crucially on a technicality 
^Jy of the definition of isotropic singularities, and not just on the uniform 

degeneracy of the metric at the singularity. 



PACS: 04.20.Dw, 98.80.H, 04.20.-q 

1 Introduction 



Isotropic singularities were introduced into the literature in |L| as a tool 
with which to study Penrose's Weyl Tensor Hypothesis and in an attempt to 
unify various studies of the initial cosmological singularity. Since then, they 
have received a lot of attention in the context of this and other hypotheses; 
in particular, the Vanishing Weyl Tensor Conjecture (that a cosmology ad- 
mitting an isotropic singularity with initially vanishing Weyl tensor evolves 
uniquely to a Robert son- Walker universe) has recently been proven for a class 
of equations of state 0. The basic idea of an isotropic singularity is that 
the singular behaviour may be removed by a conformal rescaling, where the 
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conformal factor vanishes at the singularity. The formal definition runs as 
follows |, |. 

A space-time (A4, g) is said to admit an isotropic singularity if there exists 
a space-time (Ai,g), a smooth cosmic time function t defined on M. and a 
function Q(t) such that 

(i) M. is the open submanifold t > 0. 

(ii) g ab = Q 2 (t)g ab on Xi, with g ab C 3 and non-degenerate on an open 
neighbourhood of t — 0. 

(hi) tt e C[0,b] n C 3 (0,b} for some b > 0, Q(0) = and fi(t) > for 
te(0,b}. 

(iv) The limit 

A = lim j— — / 1 

t^o+ (fi') 2 

exists, where the prime indicates differentiation with respect to t. 

The low differentiability requirement on Q at t = 0, which is identified with 
the singularity, is required to allow sufficient generality amongst this class of 
space-times. (If Q G C°° [0, b] and space-time is filled with a perfect fluid with 
a barotropic equation of state, then the fluid pressure p and density p must 
obey p = 3p (4|.) The following results may be derived using the conditions 
imposed upon Q and are essential for what follows 0: 

Lemma 1 A < 1. 



Lemma 2 lim t ^ +(^ / ) -1 ^ = 0. 



Lemma 3 (Q'^tt ~ (1 - X)t as t -> + . 

Such singularities are usually studied in a cosmo logical context, and so 
it is generally assumed that there exists a congruence of time-like curves 
in (Ai,g) representing the aggregate flow of matter in the space-time. We 
will assume this congruence to be geodesic. So let v a be the unit tangent 
to a congruence of time- like geodesies of (Ai,g), with typical member 7. 
We note that as the definition of an isotropic singularity is independent of 
the matter content of the space-time in question, this congruence need not 
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have the interpretation mentioned above. So we are studying any time-like 
geodesic (subject to one further restriction given below) which emanates from 
an isotropic singularity. We can associate with this congruence a congruence 
of time-like curves in (Ai,g) with unit tangent v a = Qv a . The congruence 
associated with v a is said to be regular if the conformal image v a is C 3 on an 
open neighbourhood of t = 0. So let (Ai,g) admit an isotropic singularity, 
and a regular time-like geodesic congruence with unit tangent v a . 

Our aim here is to study carefully the idea of the isotropy of the singu- 
larity. To this end, we show that every Jacobi field along 7 is crushed to 
zero length at the same rate as the singularity is approached. Thus physical 
objects carried along 7 are crushed to zero volume at the same rate in every 
direction; the crushing effect is isotropic. The rate at which this crushing 
takes place depends only on the index A. This is perhaps not surprising; 
indeed it verifies on a rigorous level part of what one wants to model for such 
singularities (see Tod's comments about blobs of fluid in f7|), but we shall 
see that the crushing effect arises not just from the fact that Q vanishes at 
the singularity, but that the regularity condition of the existence of the limit 
A and the assumption that t is smooth on (Ai,g a b) are crucial. As a useful 
by-product, we obtain the general asymptotic behaviour of the conformal 
factor Q as a function of proper time along 7 in the limit as the singularity 
is approached. We use the metric and curvature conventions of ||. 



2 Decomposition and solution of the geodesic 
deviation equation. 

A Jacobi field along 7 is a vector field fj a orthogonal to v a which satisfies the 
geodesic deviation equation 

D 2 f) a = -R a bcd f} c v b v d , (1) 

where R a bcd is the Riemann tensor of (Ai,g) and D = v a V a is the covariant 
derivative along 7. The solutions of ([I]) form a six- dimensional vector space 
orthogonal to v a . To solve (|l|), we use the assumptions above to identify 
the singular and non-singular terms of the right hand side. We obtain the 
asymptotic behaviour of these terms in the limit s — > + , where s is the 
proper time along an integral curve of v a , whose origin is fixed so that the 
(past) singularity is approached as s — > + . 
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By the standard formulae for conformally related metrics ||, 

R a bc<fl c v h v d = Vt- 2 R a bcd fj c v b v d 

+\(S a c Q bd - g bc Q a d - 5 a d Q bc + g bd Q a c )fj c v b v d , (2) 

where 

n ab = -4tt~ 3 V a V b tt + 8tt~ 4 V a ttV b tt - 2tt~ A V c ttV c ttg ab . 

The assumptions made above allow us to simplify the right hand side of 
(0) considerably We note first that as v a is assumed to be geodesic, then the 
relationship between the acceleration vectors of i) a and v a which comes from 
the conformal transformation 

Dv a = tt- 2 Dv a + tt- 3 Dttv a + tt- 3 V a tt 

yields 

V a tt = -DVtv a - tta a , (3) 
where D = f a V a and a b = v a V a v b . Then also 

V a ttV a tt = -(Dtt) 2 + tt 2 a b a b 

and 

7) a V a fi = - W, (4) 

and by definitions, fj a v a = fj a v a = 0. 

The scalar term in the first bracketed term of (§) can be written 

Q bd v b v d = -AQ- 1 !) 2 ^ + 2Q- 2 {DQ) 2 - 2Q- 2 a b a b , (5) 

where we have used D = VT l D when acting on scalars. 

Next, a careful calculation relying heavily on (Q) and geodesy of v a shows 
that 

Q bc fj c v b = 4Q- 2 Da b fj b . (6) 

Finally, the same conditions show that 

fl a b fj b = AQ- 2 {DQV a v b + {DQv a + a a )a b + V a a b )f) b 

+2Q- 2 {{DQ) 2 -a b a b )f) a . (7) 
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Thus from equation (^j and equations (f§)-([F|) the geodesic deviation equa- 
tion ([|) reads 

D 2 rj a + A\f) b = 0, (8) 

where 

A\ = A 1 a b + A 2 \ + A 3 a fe , 

with 

A 1 \ = Q- 2 (R a cbd v c v d - a a a b - v a Da b - V a a b ) =: Q~ 2 U^, (9) 
A 2 a b = -Q~ 2 DQ{V a v b + v a a b ) =:Q- 2 DQW b a , (10) 
A 3 a fe = -n- l D 2 Q5 a b . (11) 

We note an essential point for what follows: the definition of an isotropic 
singularity and the assumption that v a is regular, i.e. that v a is C 3 , implies 
that the tensors U b , W b are finite at the singularity. In the usual way, we can 
project (|8]) onto a parallel propagated tetrad, and so treat ff as a Euclidean 
3- vector and D as an ordinary derivative (see section 4.1 of ||). 

In order to obtain the asymptotic behaviour of solutions of (§) as s — > + , 
we must determine the asymptotic behaviour of Q and its s— derivatives (DQ 
etc.) in this limit. As a first step towards this, note that 

dq = n- l Dn = n- l n'Dt. 



By the assumption that t is a smooth cosmic time function throughout a 
neighbourhood of t = in (A4,g) and the fact that v a is C 3 , we must have 
thatQ 

t := lim Dt > 0. 
(Recall that s — > + corresponds to t — + + .) Thus 

DQ~i ^. (12) 

Similarly, the limit 

t := lim D 2 t 

s^0+ 

1 It is of course possible that Dt\t=o = 0. However the following analysis still holds in 
this case, as regardless of whether or not i$ 1 <o converges or diverges, the dominant term 
in L a b is always £3% as can be seen from equations (p"g|-p0|) in the text. 
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must exist and be finite, but may be positive, negative or zero. From these 
observations, we find that 

D 2 Q ~ Q- 2 Q'i + {Q- 2 Q" - Q- 3 Q' 2 )il. (13) 
Next, consider the term a ^ . Both numerator and denominator vanish 

' s 

in the limit s — > (see Lemma 2 above), and so we can try to use l'Hopital's 
rule to evaluate the limit of the ratio. This gives a finite and hence legitimate 
result: 

, ft 2 at , / „ n 2 n"\ « 

lim — i — = lim 2Q - — — Dt 

S -,Q+ S s^0+ \ (fi') 2 / 

/ QQ" \ 
= lim 2 - — — r Dt 

= t (2-A). (14) 

The same reasoning may be applied to the quantity lnfi/lns: 

lnfi , Dflnfi) 
lim = lim — 

s ^o+ Ins s ^o+ £)(lns) 

stt'Dt 

= hm 

where we have used flip) . This yields 

~ S^A . (16) 



Recall that A < 1, so that this power is always positive. We note that 
the relations (14|) and (pi]) hold for all isotropic singularities with a regular 
congruence of time-like curves: the assumption that 7 is geodesic was not 
used in the derivation of these equations. 

From equations (|T2l)-(|T6|) we can determine the asymptotic behaviour of 
the singular terms in (|8|). It is convenient to rewrite (§) in the form 

s 2 D 2 f} a + L a b f) b = 0, (17) 

with the obvious defintion of L a b (and its constiuents L 1 a f) etc. below). The 
singular coefficient of L\ L h is 

s 2 Q- 2 ~ s ^. (18) 
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We note that 

2 - 2A 

°< — < 2 - 

The singular coefficient of L 2 a b is 



s z Q-' 2 Dtt = -s z tt- 3 tt'Dt 



2 - A 



1 l-A 

S2=*. (19) 



Notice that this term dominates the coefficient of Li a b . The coefficient of 
L 3 a b satisifes 

-s 2 n- l D 2 n ~ -s 2 fi- 3 fi , t - S 2 (fi- 3 fi ,, -fi- 4 (fi , ) 2 )io 

2n-3n')' S (O ) / fifi \ - 2 



1 — A £n hz>: 

^s—. (20) 



(2 - A) 2 (2 - A)t 
Thus the dominant term in L a b is Ls a b , and we can write 

where M a b(s) is continuous at s = + . The s— dependence of the second term 
here indicates that L a b has a low degree of differentiability, indicating that 
the solutions of (|T?j ) cannot be written as power series in s. Nevertheless, 
the leading order behaviour of a fundamental matrix for ([17]) can be found 
by looking for a formal solution whose entries are formal logarithmic sums, 
that is functions of the form 

oo 
j,k=0 

where the fjk are formal Laurent series in s which vanish for sufficiently 
large j + k and the fij are complex numbers 0. It is easily verified that the 
leading order term in such a series corresponds to a fundamental matrix for 
the diagonal equation 

2^,2 ~a, 1 ~ A ~a q 
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A fundamental matrix for this equation is given by 



where I3 is the 3x3 identity matrix and y(s) satisfies the scalar equation 

s 2 y'\s) + -nj^ut*) = 0. (21) 

The general solution of this equation is easily found using the method of 
Frobenius. There are slight variations (the presence of logarithmic terms in 
the solution) depending on whether or not A = 2n/(n + l) for some integer 
n; however all that concerns us is that the general solution satisfies 

y(s) ~ cs u , 

where c is constant and v (positive) is the smaller of 1/(2 — A) and (1 — 
A)/ (2 — A). Thus the general solution of (|21|) satisfies y — > as s — > + . 
Hence there exists a fundamental matrix E for (|17]) which satisfies 

E{s) ~ y(s)I 3 , as s -> 0+ 

and y — > as s — > + . Therefore every solution r) a of ( [P7D shrinks to zero 
as the singularity is approached; the crushing effect of the singularity is the 
same in every direction. This reinforces the notion of isotropic behaviour at 
such singularities. We note further that the rate at which the Jacobi fields 
are crushed depends only on the index A and not on any details of the region 
of space-time from which the singularity is being approached. Thus as one 
would expect, there is an accompanying degree of homogeneity allied to the 
isotropic behaviour at the singularity. 



3 Unexpected behaviour at a singularity which 
is nearly isotropic 

To further elucidate the role played by the defining conditions of an isotropic 
singularity, we consider an example of a space-time which has an initial 
singularity with the same broad features as an isotropic singularity, but which 
does not qualify as such on a technical basis. We show that the approach to 



8 



this singularity is not accompanied by the isotropic crushing found above for 
isotropic singularities. 

We consider the spherically symmetric conformally flat space-time (so 
that g a b above is the metric tensor rjab of Minkowski space-time) with con- 
formal factor 

n{r,T)=TeKp(~Y (22) 

where n > and T, r are respectively a global Cartesian time coordinate and 
the radius function of Minkowski space-time. We take T > 0; then T = 

2 

is a past singularity of the space-time. Note that since ^ > 0, f2 — > as 
T — > + for all values of r. For T > 0, r = is the regular centre of the 
space-time, and is a time-like geodesic 7 with unit tangent 

v a = fT 1 ^. 

The equations governing the Jacobi fields along such a geodesic are given in 



[ 10[| • The two tangential Jacobi fields along the radial geodesic r = r 



constant have norms of the form 

x(s)= Cl R(s) J -^ + C2 i?( s ), (23) 

where R is the radius function of the space-time, s is proper time along 
the geodesic and are arbitrary constants giving the 2-parameter general 
solution. Along the geodesic in question, R = r Q(r ,T) and ds = QdT. 
Fixing the origin of s so that the singularity T = occurs at s = 0, we see 
that, modulo an irrelevant multiplicative constant, 



lim x(s) = lim Q(tq,T) 



dT 



s^o+ v ' t^o+ v 'J n(r ,T) 
Using l'Hopital's rule, this gives 

r £n J r \ 

lim x(s) = — lim — , = 0. 

Notice that this result applies in the case r$ = 0, and so applies to the 
geodesic 7. 

The norm y(s) of the radial Jacobi field is governed by the equation 

D 2 y(s) + Fy(s) = 0, (24) 
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where 

F(r, T) = tt- 4 (tt 2 T - Ql) - tt- 3 (tt TT - tt rr ), 
where the subscripts indicate partial derivatives. Using ds = QdT, we have 

D 2 y = n- 2 y - n- 3 n TV , 

where the overdot represents differentiation with respect to T. Then a 
straightforward calculation shows that, along 7, equation (24) may be written 



T 2 y -Ty- (2T 2 ' n - l)y = 0. 

If n < 2, the general solution of this equation vanishes at T = 0. We consider 
the case n — 2. Then the general solution may be written down explicitly; 

y (t) = Cl T 1+ ^ + c 2 T 1 -^ 2 . 

Thus the general solution, and hence the general radial Jacobi field, blows 
up at the singularity T = 0. The same holds for n > 2. So the deformation 
along this geodesic is not isotropic; infinite stretching occurs along the radial 
direction, and infinite crushing along the tangential directions. Note however 
that the characteristic volume along 7 that has norm V = \x 2 y\ satisfies 

V ~ VoT 3 -^ 2 -> 0. 

We investigate which of the conditions for an isotropic singularity this 
example fails to satisfy. We fix n = 2. Then 

* br ~ ( 2r 2 \ ^ .r 4 



rj ab Q a n b = -exp (1 + 4—). (25) 

We note therefore that V a Q is timelike for all (r, T) with T ^ 0, but is 
degenerate for Q = 0. Thus Q = does not represent a hypersurface within 
the family of hypersurfaces Q =constant. However we know that Q = <^ 
T = 0, so the singularity may be identified with the space-like hypersurface 
T = 0. (This is further verified by showing that each point (T = 0,r = r r ) 
of the hypersurface is a past endpoint for both an ingoing and an outgoing 
radial null geodesic of the space-time.) 

According to the definition of ||, the vanishing of V a f2 on Q = dis- 
qualifies this example from membership of the class of isotropic singularities. 
Similarly, it is not possible to find an appropriate time function Z = Z(Q) 
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satisfying Z(0) = and V a Z ^ 0. That is, in the language of [0, there 
is no conformal gauge choice which allows this example to be considered an 
isotropic singularity (see also [0]). This failure to qualify may be clarified 
by considering how the example fails to qualify on the basis of the definition 
above, which does not explicitly include the requirement V a f2 ^ 0. 

First, we note that taking t — O, the function t is a cosmic time function 
on Ai, i.e. t G C°(Ai) and increases along every future-directed causal curve 
of A4f\ This is clear for the open submanifolds t ^ 0, as equation (p5|) shows. 
Let s a be tangent to a future-pointing causal curve a of M.. Then s a V a t > 
for all t ^ 0, but s a V a t\t=o = 0, so certainly t is non-decreasing everywhere 
along a. But crossing t = is equivalent to crossing T = 0; this increases 
from negative to positive along a, which is future directed. Hence t also 
increases from negative to positive along a. So t does indeed increase, in the 
strict sense, along every future-directed causal curve of M.; £|<t(t 2 ) > ^o-(n) 
for values T2 > T\ of a parameter r along cr which increases into the future. 
In particular, we have Dt = 1 along the geodesic r = 0, and so t = 1 for 
this example. 

With t — Q, conditions (i)-(iv) above are trivially satisfied, with A = 0. 
Also, the geodesic 7 is regular with geodesic conformal image. Thus it is 
only the failure of t := Q(r,T) to be a smooth cosmic time function on 
(Ai,r] ab ) which precludes this example from satisfying all the hypotheses 
of an isotropic singularity. In fact t fails to be smooth (C°°) only at the 
point (r, T) = (0,0). Therefore this example shows the importance of this 
assumption for isotropic singularities. In terms of the relationship with the 
definitions of f| and @, as pointed out by Tod in 0, the current definition 
CO Ml) assumes implicitly that among the allowed cosmic time functions of 
the definition above, one - t say - can be chosen so that V a t 7^ at the 
singularity. Assuming smoothness for a particular choice t would allow one 
to find the "correct" t = i(t) (this issue may deserve closer attention). The 
lack of smoothness of our t in this example indicates that no such i exists; 
the degeneracy of Q at the singularity, i.e. V a Q = is gauge-invariant. 



2 This is the definition of |T^]. Note the inequivalence to the definition that V a t is time- 
like everywhere. The former definition is more general, allowing (i) less smooth functions; 
(ii) examples such as t = T 3 , where T is a global Euclidean time co-ordinate on Minkowski 
space-time. 
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4 Conclusions 



We have studied some technicalities of the definition of isotropic singularities. 
As one would expect for such a definition, the behaviour at the singularity 
is indeed isotropic, in the sense that Jacobi fields carried along a geodesic 
running into such a singularity are crushed uniformly in every direction. 
This verifies in a rigorous manner what one intends to model with such a 
definition; a blob of matter would be crushed uniformly to zero volume as the 
singularity is approached. We have also found homogeneous behaviour at the 
singularity; the rate at which the Jacobi fields are crushed depends only on 
the index A of the definition above, and not on the region of space-time from 
which the singularity is approached. Thus although space-times containing 
isotropic singularities are usually inhomogeneous |J, the inhomogeneity is 
washed out in the approach to the singularity. 

We have derived this result for time-like geodesies subject to a certain reg- 
ularity condition. For non-geodesic motions, the deviation equation becomes 
more complicated (see equation (4.8) of ||). However the generalisation of 
the result to bounded acceleration curves (for example the fluid flow lines 
in the examples of Jl2j ) does not appear to be inconceivably difficult: since 
in this case equation (|3p remains true in the limit Q — > 0, the asymptotic 
analysis should follow through. 

The example of section 3 above shows how these conclusions depend cru- 
cially on a technicality of the definition of |3j, namely that the conformal 
factor is a function of a cosmic time function which is smooth on the back- 
ground space-time. In our example, the cosmic time function fails to be even 
C 1 throughout the hypersurface t = 0. This manifests itself in the fact that 
regardless of how one tries to choose the conformal factor Q and appropri- 
ate time function Z (see [0]) for this particular space-time, one finds the 
degeneracy V a Z = at the singularity. 
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